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ABSTRACT 


This  report  presents  an  analysis  of  the  axisymmetric  elastic  and 
plastic  stresses  and  deformations  in  thick  wall  cylindrical  shells  sub¬ 
jected  to  internal  dynamic  pressures.  A  direct  numerical  approach 
called  the  discontinuous -step  analysis  is  used.  This  analysis  is  based 
on  the  direct  use  of  the  boundary  conditions  and  the  applicable  physical 
laws  to  propagate  dynamic  changes  in  the  cylinder  by  finite  steps.  Re¬ 
flection  of  stress  waves  from  both  inner  and  outer  boundaries  is  auto¬ 
matically  generated.  The  validity  of  the  method  is  checked  by  compar¬ 
ison  of.  numerical  results  in  the  elastic  range  with  published  results 
for  thick  wall  cylinders.  Comparison  is  made  with  experimentally 
measured  strains  from  the  high-pressure  section  of  a  hypervelocity 
launcher.  This  analysis  assumes  that  the  work  hardening  of  the  mate¬ 
rial  is  independent  of  the  strain  rate  and  is  constant  for  a  large  varia¬ 
tion  of  plastic  strain.  Stress -strain  relationships  are  derived  for  the 
condition  of  plane  stress  in  the  cylinder  which  is  held  to  be  represent¬ 
ative  of  the  actual  conditions  in  the  launcher  high-pressure  section. 

The  digital  computer  program  developed  from  this  study  predicts  the 
distribution  of  dynamic  stress  and  strain  throughout  the  cylinder,  the 
internal  radial  growth,  the  distribution  of  particle  displacement,  the 
distribution  of  yield  stress  in  an  autofrettaged  cylinder,  and  the  resid¬ 
ual  stress. 


iii 


AEDC-TR-70-42 


CONTENTS 

Page 

ABSTRACT . iii 

NOMENCLATURE . vi 

I.  INTRODUCTION .  1 

II.  BASIC  EQUATIONS 

2.  1  Elastic  Range . 1 

2.  2  Plastic  Range . 2 

2. 3  Momentum  Balance .  4 

2.4  Stress  Propagation .  5 

III.  NUMERICAL  PROCEDURE .  6 

IV.  RESULTS  AND  DISCUSSION .  10 

V.  RECOMMENDATIONS .  12 

REFERENCES .  12 

APPENDIXES 

I.  ILLUSTRATIONS 
Figure 

1.  Nomenclature  of  Cylindrical  Cross  Sections . 17 

2.  Stress -Strain  Curve  for  Ductile  Steel . 18 

3.  Force  Balance  on  Semicircular  Element . 19 

4.  Response  of  Cylinder  to  Step  Input  PQ.  Comparison  of 

Solution  by  Characteristics  (Chou  and  Greif)  and  Discon¬ 
tinuous  Step  Method . 20 

5.  Pressure  versus  Time  Distribution  Approximating  to  Shot 

T  1237  (221,  000-psia  Peak  Pressure) . 21 

6.  Pressure  versus  Time  Distribution  Approximating  to  Shot 

T  1247  (194,  000-psia  Peak  Pressure) .  22 

7.  Comparison  of  Measured  and  Calculated  Strain  for  Shot 

T-1237  .  23 

8.  Comparison  of  Measured  and  Calculated  Strain  for  Shot 

T-1247  .  24 

9.  Predicted  Radial  Distribution  of  Tangential  and  Radial 
Residual  Stress  after  Lower  Pressure  Shot  (194,000  psi) 

and  Higher  Pressure  Shot  (221,000  psi) .  25 


v 


AEDC-TR-70-42 


Page 


II.  TABLES 

I.  Radial  Growth  (DELR) .  26 

II.  Peak  Tangential  Strain .  26 


NOMENCLATURE 


Constant 

Area,  ith  element  sq  in. 

Inner  radius,  in. 

Outer  radius,  in. 

Elastic  wave  velocity,  in.  /  sec 
Wave  velocity  elastic  on  plastic 
Plastic  wave  velocity,  in.  /  sec 
Elastic  modulus,  psi 
Plastic  modulus,  psi 
Force,  lb 

Gravitational  constant,  in. /sec2 
Bulk  modulus,  psi 
Element  length,  in. 

Mass,  lb 

Transient, internal  pressure,  psi 
Internal  pressure,  psi 
Step  input  pressure 
Radius,  in. 

Time,  sec 

Radial  growth  ith  element,  in. 

Particle  displacement,  in. 

Particle  velocity,  in.  / sec 

Total  time  interval  of  calculation,  sec 
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Y  Original  yield  point,  psi 

6  Arbitrary  small  constant,  psi 

e  Strain,  in.  /  in. 

v  Poissons  ratio 

O 

p  Density,  lb/ in. 

cr  Stress,  psi 

or0  Work  hardened  yield  point,  psi 

t  Shear  stress,  psi 


SUBSCRIPTS  AND  SUPERSCRIPTS 

e  Elastic  range 

i  Index  for  identification  of  radial  grid  coordinate 

m  Total  number  grid  coordinates 

p  Plastic  range 

r  Radial  direction 

x  Arbitrary  radial  direction 

z  Longitudinal  direction 

0  Tangential  direction 
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Mehta  and  Davids  (Ref.  1)  have  investigated  the  propagation  of  tran¬ 
sient  stress  waves  in  cylinders  subjected  to  various  types  of  dynamic 
loading.  The  study  was  limited  to  the  elastic  range  of  stress,  and  the 
authors  employed  their  discontinuous -step  analysis.  A  force  balance 
and  the  stress-strain  relations  were  formed  as  the  basic  governing  phe¬ 
nomenon.  These  were  used  together  with  the  proper  boundary  condi¬ 
tions  to  predict  discrete  changes  in  stress  in  finite  time  steps.  With 
slight  alteration,  this  is  the  method  employed  in  the  present  study. 

There  are  significant  differences  in  the  analysis  of  dynamic  stress 
propagation  for  the  plastic  range  as  opposed  to  the  elastic  range.  In  the 
elastic  state,  the  material  retains  constant  properties  of  modulus,  and 
the  characteristic  is  well  defined.  In  the  plastic  range,  the  modulus 
varies  with  the  state  of  stress,  and  wave  velocity  changes  from  that  of 
the  elastic  state.  With  the  problem  at  hand,  that  of  predicting  the  tran¬ 
sient  stress  in  the  cylindrical  high-pressure  section  of  a  hypervelocity 
launcher,  both  the  elastic  and  plastic  states  of  the  material  jointly  exist 
in  the  higher  pressure  regimes.  Any  stress  analysis  must  give  cogni¬ 
zance  to  the  fact  that  the  cross  section  of  the  cylinder  first  becomes 
plastic  at  the  core  and,  as  pressure  increases,  the  plastic  boundary 
advances  radially  outward. 

Davids,  Mehta,  and  Johnson  (Ref.  2)  have  investigated  the  propaga¬ 
tion  of  transient  stress  waves  in  the  el asto -plastic  regime  for  a  spherical 
body.  In  this  case  also,  the  authors  applied  their  discontinuous -step 
analysis  with  emphasis  on  the  solution  of  the  radial  displacement  at  the 
outside  boundary  of  the  sphere.  An  experiment  was  performed  by  the 
detonation  of  a  charge  placed  within  the  cavity  of  a  thick-wall  aluminum 
sphere.  Analytical  and  experimental  results  did  not  agree  well,  and  the 
authors  concluded  that  the  analysis  was  only  applicable  for  cases  of  load¬ 
ing  of  a  nondecreasing  nature  (i.  e. ,  the  rising  portion  of  a  pressure 
pulse).  As  will  be  indicated  later,  the  present  study  removes  this  limi¬ 
tation  from  the  discontinuous  method. 

SECTION  II 
BASIC  EQUATIONS 


2.1  ELASTIC  RANGE 

Consider  a  cross  section  of  the  cylindrical  body  as  shown  in  Fig.  la 
(Appendix  I).  The  radial,  tangential,  and  longitudinal  stresses  (crr,  erg, 
and  ctz)  load  an  element  as  shown  in  Fig.  lb.  In  the  problem  at  hand,  the 
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internal  impulsive  pressure  loading  is  caused  by  a  piston  moving  along 
the  longitudinal  axis  of  the  cylinder  and  compressing  gas  trapped  therein. 
The  principal  loading  is  then  radial  in  nature,  and  longitudinal  strain  is 
unimportant  to  the  analysis.  It  is  held  that  conditions  of  plane  strain 
exist;  hence,  ez  is  zero.  In  the  elastic  range  the  stress-strain  relations 
are  then 


iq  <=  sr  [oq  -  v{oz  +  £Tr)] 

L  1 

«r  =  -  v(oq  +  az)] 

(z  =  0  =  —  [<TZ  -  v\.OQ  +  (Jr)] 

El 

Substituting  for  crz  from  Eq.  (3)  yields 

£  J 

[(1  “  +  ^ 


(1) 

(2) 

(3) 

(4) 


and 


°Q 


(l  +  i/)(l-2i/)  +  £r 


(5) 


Equations  (4)  and  (5)  relate  the  radial  and  tangential  stress  to  the 
strains. 


2.2  PLASTIC  RANGE 

Following  Hopkins  (Ref.  3,  page  128),  a  work. hardening  model  is 
assumed  for  the  stress-strain  relation  of  the  material.  Let  the  stress- 
strain  curve  of  the  steel  in  uni- axial  loading  be  as  shown  in  Fig.  2. 

Here  the  effects  of  rate  of  strain  and  rate  of  hardening  are  neglected. 
The  plastic  modulus  E2  is  taken  as  constant  in  this  assumption  of 
linear  strain  hardening. 

The  plastic  strains  for  the  small  cylindrical  segment  are  taken  to 
be  incompressible  (see  Hopkins,  Ref.  3,  page  113).  Then,  for  plain 
strain, 

er  ^  $  =  0  (6) 

but  the  total  strain  is  the  sum  of  the  plastic  and  elastic  strain, 

fr  =  er  +  4r  and  e6  =  +  e§  (7) 
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Combining  Eqs.  (6)  and  (7),  we  have 

(r  -  <6  =  eT  +  e6  (8) 

Now,  from  Eqs.  (1)  and  (2)  for  the  plane  strain  case 

e*  =  [ar(l  -  i/)  -  vo$) 

M 

<CQ  =  -  v)  -  VOr] 

Substituting  Eqs.  (9)  and  (10)  into  (8),  one  finds  the  bulk  compress¬ 
ibility  relation  for  the  cylindrical  plain  strain  case  as 

oT  +  oe  =  +  £0')  (11) 

Now  for  the  elastic  state,  one  may  assume  a  simple  Tresca  yield 
criterion  (point  A,  Fig.  2).  The  principal  strain,  depicted  in  Fig.  3a, 
exists  in  the  r-9  plane  and  the  yield  criterion  is 

n  ~  *r  =  Y  (12) 

Realizing  that  this  criterion  attains  only  in  the  elastic  range  or 
for  a  perfectly  plastic  material,  one  sets  to  the  task  of  modifying  the 
criterion  to  relate  stress  and  strain  in  the  plastic  range  (point  C, 

Fig.  2).  The  yield  criterion  is  of  the  form 

a  =  Y '  +  E2  (f)  (12a) 


(9) 

(10) 


where  a  is  the  principal  stress  and  e  is  the  equivalent  total  strain 
(Ref.  4).  Under  these  assumptions,  Eq.  (12a)  becomes 


OQ  -  O' 


(12b) 


Combining  Eqs.  (11)  and  (12b),  one  obtains  expressions  for  stress  in 
terms  of  total  strain  in  the  plastic  state.  Thus, 


and 


Ox 


(3K-E2)  (3K  +  E2) 

- iQ  +  - 

2{l  +  i/)  2{l  +  i/) 


(13) 


_  J3K +_£,)  (3K  -  E2) 

2(1  +i/)  2  (1  +  i/) 


Y/2 
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2.3  MOMENTUM  BALANCE  o 

The  essence  of  the  discontinuous -step  method  consists  of  defining 
the  force  difference  on  a  given  element  from  the  known  instantaneous 
stress  distribution.  In  a  given  time  increment,  this  force  creates  an 
impulse  (force  x  time)  giving  rise  to  a  change  in  momentum  of  the 
element  (mass  x  change  in  velocity).  Considering  the  semicylindrical 
element  of  Fig.  3,  the  net  force  acting  in  a  given  direction  is  (Ref.  1) 

Fx  =  [(ari+1  (2r,+  l0  -  (2rj £)  -  ^  (2r-,+ i  -  2rj)£] 

Now,  the  inner  surface  area  of  the  ith  element  is  given  by 

Ai  =  ffTi? 


and  the  net  force  may  be  written 

Fx  =  4  tw+1  Ai+l  -  ct;  Ai)  -  O0i( Ai-!  -  Ai)] 

Now,  the  impulse -momentum  relations  equate  the  total  impulse  Fxdt 
to  the  momentum  m^dv  of  the  mass  center.  The  mass  center  of  the 
semicylindrical  element  of  radius  R  is  located  on  the  x  axis  (Fig.  3) 
at  R  -  2B/w.  Then 


Mjdv  =  Fxdt 


and 


dv  2  dv 
dt  77  dt 

Substituting  in  the  expression  for  Fx  and  canceling  2/tt,  we  have 

dv  =  [(<7r,+  1  Aj+i  -  ffr'Ai)  +  (70  (A;  -  Ai+1)]dt/M.  (15) 

where  Mj,  the  mass  of  the  ith  element  is 


By  use  of  Eq.  (15),  one  can  define  the  incremental  change  in  particle 
velocity  dv  that  occurs  in  a  given  time  increment  as  the  result  of  the 
radial  and  tangential  stress  distribution. 
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2.4  STRESS  PROPAGATION 

As  indicated  in  Fig.  3,  one  first  divides  the  cylinder  into  a  finite 
number  of  elements  im.  This  radial  grid  is  then  fixed  in  time,  and 
the  cells  within  the  grid  assume  instantaneous  states  of  strain  in  the 
elastic  or  plastic  region  depending  on  the  internal  pressure  and  the 
elapsed  time.  The  state  of  any  cell  can  be  defined  at  any  time  by 
assigning  the  variables  of  particle  displacement  u,  particle  velocity  v, 
strain  e  ,  and  stress  a.  Additionally,  the  state  of  any  cell,  of  width  dr, 
can  be  changed  over  a  time  increment  dt  by  defining  the  incremental 
change  in  displacement  du,  velocity  dv,  and  strain  de  . 

In  the  radial  coordinate  set  used,  the  tangential  and  radial  strains 
are  by  definition  (Ref.  5), 

u  ,  du 

*6  =  T  and  fr  =  7T 


In  starting  the  iteration  process,  an  impulse  of  pressure  P(t)  is  given 
to  the  core  of  the  cylinder  in  the  first  time  increment.  This  gives 
rise  to  a  stress  wave  having  ther  elastic  wave  velocity 

/  E  (l-*4  \* 

Le  =  \J  (i-vUT-w)  (16) 

or  the  plastic  wave  velocity  (Ref.  3,  Eq.  6.40) 

„  /K(l  +  E2/3K>Y* 

P  -  E2/9K )p )  (17) 

depending  on  the  magnitude  of  the  principal  stress. 

In  a  given  time  t,  the  state  of  the  system  is  examined  at  each  indi¬ 
vidual  cell.  If  the  radial  velocities  of  two  adjacent  cells  i  and  i  +  1 
are  v^  and  v^  +  respectively,  from  the  basic  definition  of  radial  strain 
one  can  define  an  average  radial  strain  increment  de  during  the  time 
interval  dt  by 

d^r  =  (vi+  1  ~  vi)  >  *  =  1,  2,  .  •  ,  im  (18) 

where  dr  is  the  constant  radial  grid  increment.  Since  the  element  may 
already  be  strained  at  the  beginning  of  the  time  interval  dt  by  an  amount 

eri  +  \  the  new  strain  level  at  the  end  of  the  increment  i+*^'  at  the 

end  of  the  interval  is  the  cumulation 

(fri+  I)'=  cri+  1  df  r 


5 


AEDC-TR-70  -42 


Then,  the  state  of  strain  at  the  end  of  the  time  interval  is  reinitiated  by 
replacing  the  old  level  by  the  new  level,  or  using  the  replacement 
symbol, 


+  *ri+l  +  dfr 


(19) 


Similarly,  from  the  mathematical  definition  of  tangential  strain, 

,  dt 

**-Vi  T 


(20) 


and  the  corresponding  cumulation  condition  is 

(id*'b  +  ied  (21) 

Finally,  the  incremental  change  of  velocity  dv  occurring  over  a  time 
interval  is  defined  from  Eq.  (15),  and  velocity  is  cumulatively  updated 
with 


vi  *■  vj  +  dv 


(22) 


Relations  (19),  (21),  and  (22)  form  the  essential  basis  for  the 
discontinuous -step  analysis.  One  needs  only  to  direct  these  operations 
repetitively  for  i  =  1  to  i  =  im  in  a  given  time  interval,  then  update  the 
time 


t-  t  +  dt 

continuing  over  the  time  span  of  interest.  Additional  statements  are 
necessary  to  adjust  the  strain  increments  when  the  plastic  range  is 
reached  and  to  control  the  selection  of  the  constant  radial  mesh  width  dr. 
One  must  also  be  aware  of  the  sequential  ordering  of  operations  which  is 
of  critical  importance.  These  final  contingencies  will  now  be  discussed. 


SECTION  III 

NUMERICAL  PROCEDURE 


The  first  choice  facing  the  analyst  is  in  the  selection  of  the  grid 
mesh  width  dr  and  the  time  increment  dt.  Contrary  to  the  practice 
adopted  by  Mehta  and  Davids  (Refs.  1  and  2)  the  present  study  fixes  the 
time  increment  dt  and  the  mesh  width  dr  from  the  outset  and  retains 
these  two  dimensions  as  fixed  regardless  of  the  state  of  stress,  either 
plastic  or  elastic.  The  usual  characteristic  assumption  is 

dr  =  C  jdt 
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where  Cj  is  the  appropriate  wave  velocity  for  the  elastic  or  plastic 
state.  It  should  be  noted  that  plastic  wave  velocity  is  less  than  elastic. 
For  a  fixed  time  interval,  the  radial  increment  dr  would  change  (drop) 
when  the  material  became  plastic,  leading  to  confusion  in  previously 
accumulated  stress  based  on  the  original  grid  mesh. 

Previous  investigations  (Refs.  6  and  7)  concerned  with  wave 
phenomena  involving  hyperbolic  equations  have  shown  that  computations 
are  stable  if 


dr  >  Cjdt 

From  this  fact,  an  appropriate  time  increment  is  chosen,  and  a  radial 
increment  dr  is  set  such, that  it  is  greater  than  the  elastic  wave  (fastest) 
advance  in  this  preselected  time  interval.  For  computational  purposes, 
we  have 


dr  =  A  Cedt 


where  A  is  a  constant  greater  than  unity,  such  that  dr  >  (dr)e  >  (dr)p. 

The  computer  program,  written  for  this  study,  has  the  following 
sequence  of  operations: 


(1)  Specify  given  data:  aQ,  bQ,  f,  p,  Ej,  E2,  v,  dt, 

(2)  Specify  radial  distribution  of  yield  point,  cr0  and 
initial  values  of  or  and  <jq,  er  and  eg,  if  auto- 
frettaging  exists. 

(3)  Define 


dr  =  A  Cedt 


ri  +  1  =  ri  +  dr 
A;  =  jrrj£ 

Ms  =  i  - 

(4)  Define  pulse  input  P  =  P(t) 

then  cr*  =  -  P  (boundary  condition) 

(5)  Propagation  procedure:  repeat  steps  A  through  L  for 
i  =  1  to  i  =  im. 
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(vi+i  -  Vi)dt 

A.  Increment  radial  strain  der  =  — 


dr 


B.  Accumulate  radial  strain 


-»  frl+1  +  d«r 


C.  Define  radial  and  tangential  stresses 
(1)  Elastic  State 


i+l 


(1 


- - — — 2\  kr'  +  1  (  1  -  v)  +  Vf’j 

-  v  -  2ir)  & 


°e 


(1  -  v-  2v2) 


[vfri+1  +  (i  -  v)  «y 


(2) 


Plastic  State 


ri+  1  = 


(3K  -  E2) 
2(1  +  v) 


f0  + 


GQ  = 


(3K  +  E  2 )  i 

-  tn  + 

2(1  +  u)  0 


(3K+E2)  .+  1 

-  £if 

2(1  +  v) 


(3K  -  E2) 
2(1  +  v) 


+ 


D. 


Apply  boundary  corrections 1 :  for  i  =  1  and  i  =  im  only 
(1)  Elastic  State 


o'q  =  — y  tE  i  f +  vd  +  v)oT'] 

(2)  Plastic  State2 


G0  = 


1Required  because  of  boundary  condition  4  and  defined  from  Hooke's 
Law,  Eq.  (1). 

o 

^Defined  from  bulk  compressibility  relation  and  plastic  yield  crite¬ 
rion,  Eqs.  (11)  and  (12b). 
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E.  Apply  boundary  correction:  for  i  =  im  only 

i  +  l  p. 

crr  =0 

F.  Define  tangential  strain  increment 


v:dt 

-T- 


G.  Cumulate  tangential  strain 


fr  d£(? 


H.  Define  velocity  increment 


dv  =  [(crj+1  Ai  j.  1  -  o[  Ai)  +  <70  (A;  -  Ai+l)l  — 

Mi 


I  .  Cumulate  particle  velocity 

Vi  vj  +  dv 

J  .  Define  radial  displacement  increment 

dur  =  vjdt 

K.  Cumulate  radial  displacement 


U1  «-  U*  +  dur 
r  r 


L.  Increment  Time 


l«-  t  +  dt 

Steps  A  through  L  are  repeated  for  a  specified  time  interval  to 
t  =  Xm.  If  the  induced  longitudinal  stress  is  desired  for  this  plane 
strain  case,  it  can  be  found  from 


<7‘  =  k(<71  +  1  +<7]j) 

z  r  v 

A  particular  procedure  is  followed  when  stress  levels  reach  the 
plastic  state  as  defined  by  the  usual  Tresca  yield  criterion  (see  Ref.  2, 
page  131).  Referring  to  Fig.  2,  the  state  of  stress  represented  by 
point  B  is  not  possible  physically.  When  a  new  strain  increment  de  r 
at  some  time  interval  indicates  that  the  plastic  condition  attains,  stress 
is  now  calculated  by  step  5C(2)  and  the  stress -strain  relationship  is 
established  by  the  line  AC. 
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At  this  point  in  the  analysis,  judgment  is  made  on  the  amount  of 
work  hardening  which  has  attained.  The  total  principal  stress  is 
investigated,  and  the  yield  point  Y  is  updated  to  a  new  yield  point  ffQ 
defined  by 

i  /  i  _  i  +  lv 

°o  =  iog  -  ) 


Thus,  a  radial  array  of  work  hardening  is  one  of  the  results  of  a  given 
computer  run.  As  strain  is  relieved  with  decreasing  internal  pressure, 
elastic  relief  is  allowed  from  point  C,  Fig.  2.  Additionally,  it  is 
possible  to  prestore  a  radial  distribution  of  work  hardened  yield  points 
cri  and  residual  stresses  and  strains  at  the  outset  of  computations. 

This  is  representative  of  an  autofrettaged  cylinder  (step  2).  The  final 
distribution  of  residual  stress  and  strain,  work  hardened  yield  points 
and  displacements  can  be  retained  for  subsequent  transient  pressure 
loadings  of  the  cylindrical  section  if  required. 


SECTION  IV 

RESULTS  AND  DISCUSSION 


In  order  to  check  the  accuracy  of  the  derivations  and  the  program, 
comparison  was  made  with  published  solutions  for  thick  wall  cylinders. 
The  literature  only  contains  solutions  limited  to  the  elastic  range;  hence, 
this  part  of  the  analysis  could  be  verified. 

Figure  4  gives  the  solution  for  dimensionless  tangential  stress  at  the 
inner  radius  aQ  for  a  cylinder  of  outer  radius  bQ  =  2a0.  The  core  is 
given  a  step  input  of  pressure  at  time  zero.  The  solution  from  this  study, 
using  the  discontinuous -step  analysis,  is  compared  with  that  of  Chou  and 
Greif  (Ref.  7)  obtained  from  a  method  of  characteristics.  The  interested 
reader  should  refer  to  Fig.  4,  page  1071  of  Ref.  7  in  which  the  authors 
compare  the  results  of  Mehta  and  Davids  (Ref.  1)  with  their  solution. 

The  present  study,  using  the  same  discontinuous -step  analysis,  shows 
much  improved  results  over  those  of  Mehta  and  Davids  (Ref.  1)  shown  in 
Fig.  4  of  Ref.  7.  Correction  of  an  error  in  Eq.  16e  of  Ref.  1  and  revised 
rules  on  selection  of  mesh  width  dr  are  held  to  account  for  the  improve¬ 
ment. 
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In  order  to  load  the  cylinder  internally,  it  is  necessary  to  know  the 
transient  gas  pressure  in  the  high-pressure  section  of  the  launcher. 
Pressures  here  are  so  severe  as  to  preclude  measurement  with  conven¬ 
tional  transducers.  Failing  to  know  the  pressure-time  signature  from 
experiment  for  any  given  launch  cycle,  one  uses  predicted  pressures. 
These  predicted  pressures  are  obtained  from  a  real  gas  analysis  of  the 
cycle  after  DeWitt  and  Cable  (Ref.  8).  Figures  5  and  6  show  plots  of 
pressure  from  analysis  which  corresponds  best  to  measured  piston 
velocity  and  launch  velocity  for  two  shots  of  a  hypervelocity  launcher  at 
AEDC.  Each  transient  pressure  load  was  approximated  by  a  series  of 
straight  lines  showmin  the  figures.  This  allowed  the  determination  of 
instantaneous  pressure  as  a  function  of  time  from  a  subroutine  algorithm. 

Figures  7  and  8  compare  the  measured  tangential  strain,  e  q,  at  the 
periphery  of  the  high-pressure  section  of  the  launcher  with  the  calculated 
strain.  The  strain  gages  were  located  in  the  region  of  highest  loading. 
The  transient  internal  pressures  were  taken  from  Figs.  5  and  6,  respec¬ 
tively.  Comparison  shows  that  predicted  peak  strain  agrees  well  with 
peak  measured  strain.  The  time  duration  of  the  pressure  pulse  in  the 
analysis  corresponds  with  that  of  the  measured  transient  strain.  In 
each  case,  the  computed  results  for  the  autofrettaged  cylinder  are  in 
better  agreement  with  measured  results.  The  reader  should  consider 
that  measurements  were  made  on  a  fully  autofrettaged  section. 

Figures  7  and  8  show  a  high  frequency  ripple  superimposed  on  the 
general  transient  response.  Investigation  shows  that  this  frequency 
corresponds  very  well  with  the  fundamental  ring  or  hoop  frequency  for 
the  dimensions  of  the  cylinder.  Its  period  is  significantly  that  time 
required  for  a  transient  pulse  to  be  transmitted  to  the  outer  boundary 
and  reflected  to  the  core. 

A  study  was  made  of  the  effects  of  autofrettaging  on  the  thick-walled 
radial  section.  Figure  9  compares  the  computed  residual  stress  remain¬ 
ing  in  the  hypothetical  study  of  two  launch  cycles  of  increasing  severity. 
The  subtlety  here  is  that  the  lower  pressure  cycle  (194  x  10^  psi)  is  first 
applied  to  a  virgin  section.  Wall  yielding  occurs,  the  transient  pressure 
disturbance  subsides,  and  the  new  equilibrium  state  of  static  radial  and 
tangential  residual  stress  (solid  lines)  is  predicted  from  the  basic  force 
balance.  This  radial  stress  distribution  is  prestored  prior  to  the  sub¬ 
sequent  computer  analysis  of  a  higher  pressure  launch  cycle  (221  x  10^ 
psi).  Figure  9  shows  the  new  computed  residual  stress  remaining  after 
the  second  launch  cycle.  As  expected,  the  cylinder  is  further  auto¬ 
frettaged  and  peak  radial  and  tangential  residual  stress  are  increased. 

The  first  launch  cycle  yielded  the  wall  through  42  percent  of  its  thickness. 
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The  second  shot  increased  the  plastic  radius  to  52  percent  of  the  wall 
thickness.  This  corresponds  well  with  the  location  of  the  abrupt  slope 
change  in  the  ffg  curves  of  Fig,  9. 

Table  I  (Appendix  II)  lists  the  predicted  bore  growth  of  the  cylinder 
from  the  analytical  pressure  pulses  of  shots  T-1237  and  T-1247  (Figs.  5 
and  6,  respectively).  The  table  compares  the  growth  for  a  virgin  section 
and  for  an  autofrettaged  section.  In  each  case,  the  autofrettaged  section 
showed  less  radial  growth.  Table  II  compares  the  peak  external  tangen¬ 
tial  strains  for  each  of  the  shots  listed.  It  is  seen  that  previous  strain 
history  directly  affects  the  dynamic  response.  The  autofrettaged  section 
showed  significantly  less  peak  strain.  It  is  held  that  the  large  radial 
region,  having  higher  yield  points  from  work  hardening,  remains  elastic 
for  higher  pressure  excursions.  Hence,  the  autofrettaged  section  sees 
less  strain  induced  for  the  same  internal  pressure  distribution.  The 
experimental  measurements  are  in  good  agreement  with  the  prediction 
for  the  autofrettaged  section. 


SECTION  V 
RECOMMENDATIONS 


In  order  to  provide  for  the  effects  of  strain  rate  on  the  stress -strain 
characteristics  of  the  steel,  the  present  mathematical  model  should  be 
expanded.  A  Voigt  model,  or  perhaps  a  combinational  Zener  model, 
could  be  considered  in  the  momentum  balance  of  the  radial  element. 

Thus 


M  dv  =  F  dt  -  tj£  dt 

where  rj  is  the  appropriate  loss  coefficient  for  the  material. 
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Fig.  2  Stress-Strain  Curve  for  Ductile  Steel 
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Fig.  4  Response  of  Cylinder  to  Step  Input  PQ.  Comparison  of  Solution  by  Characteristics 
(Chou  and  Grief)  and  Discontinuous  Step  Method 
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Fig.  5  Pressure  versus  Time  Distribution  Approximoting  to  Shot  T  1237  (221,000-psia  Peak  Pressure) 


ig.  6  Pressure  versus  Time  Distribution  Approximating  to  Shot  T  1247  (194,000-psia  Peak  Pressure) 
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Fig.  7  Comparison  of  Measured  and  Calculated  Strain  for  Shot  T -1237 
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TABLE  \ 

RADIAL  GROWTH  (DELR) 


Shot  Number 

T-1247 

T - 1237 

Non-Autofrettaged  5.5  x  10“3  in. 

9  x  10"3  in. 

Autofrettaged 

0 

TABLE  11 

PEAK  TANGENTIAL  STRAIN 

0 

Shot  Number 

T-1247 

T-1237 

Non-Autofrettaged 

1600  x  10"®  in.  / in. 

2098  x  10'6  in.  /in. 

Autofrettaged 

985  x  10~®  in.  / in. 

1179  x  lO-^  jn>  /in> 

Experimental 

1030  x  10"®  in.  / in. 

1020  x  10"6  in.  /in. 
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